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1. Introduction. 



The bivariate delayed Poisson process is a particular type of bivariate 
point process, that is a point process with two types of event occurring along 
a time axis. Bivariate point processes were first studied by Cox and Lewis 
(1972) in their Berkeley Symposium paper and here their analysis of the 
bivariate Poisson processes (without Poisson noise) is carried further. This 
bivariate point process is formed by subjecting the events of a main Poisson 
process to independent pairs of delays; the two types of delay for each main 
event produce the two types of event in the bivariate point process. Since it 
is well known that a Poisson process whose events are independently and identi- 
cally displaced remains a Poisson process, the events of each type considered 
separately form two Poisson processes. They cannot, however, be independent 
Poisson processes because their events are associated in pairs by the displace- 
ment centres (the events of the original Poisson process). The resulting 
bivariate point process is by definition a particular type of bivariate 
Poisson process because its marginal processes for events of each type are 
themselves (univariate) Poisson processes. Cox and Lewis (1972) assumed a 
model in which the delays within pairs were independent and non-negative, and 
discussed the transient and asynchronous bivariate counting processes. We shall 
mainly assume that the delays within pairs are dependent and two-sided and will 
first obtain the stationaly initial conditions. These lead to the joint distri- 
bution of the numbers of delays of each type in progress at an arbitrary time 
and a general form of the asynchronous bivariate counting distribution. Finally, 
we obtain the semi-synchronous marginal interval distribution, that is the 
distribution of the time from an event of one type to the next event of the 
opposite type; such distributions are of obvious use in the statistical analysis 



of bivariate point processes, and have been studied in the context of the 
general theory by Wisniewski (1972). Milne (unpublished Ph.D. thesis), has 
shown that the bivariate Poisson process studied here is the most general form 
for an infinitely divisible bivariate point process. 

2. Notation and Description of the Model. 

The bivariate delayed Poisson process is constructed from a main Poisson 

process of rate y. Associated with each of its events is a pair of delays 

and Y^; thus an event at time t in the main Poisson process produces 

an event at time ^ a-process and an event at time t 4- Y^ in 

the b-process . We will assume that the random variables Y^ and Y^ with 

ranges (-",") have finite means and a dependent joint distribution without 

any atoms. The pairs (Y^,Y^) associated with different main events are 

assumed to be independently and identically distributed. 

The random variables N^(t,t+t^), Nj^(t,t‘ft^) denote the numbers of 

type-a and type-b events in the intervals (t,t+t ] and (t,t+t, ], respectively 

a 0 

when there can be no confusion these will be abbreviated to 

with joint probability generating function (j .p.g. f . ) ,cp (z^, Z 2 ; t^, t^) • 

3. Stationary Initial Conditions. 

The purpose and methodology of stationary initial conditions for uni- 
variate point processes have been described in Lawrance (1972); they carry 
over quite obviously to bivariate processes and so will not be laboured here. 

We wish to study the behaviour of the process at an arbitrary time, and this 
is defined as the time t as t It is apparent for the bivariate delayed 

Poisson process that at any time t there may be a (random) number of events 
of each type in the course of being delayed. The initial conditions must 
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therefore specify jointly the distribution of these numbers and the individual 
residual delay times measured from t. We begin by dealing with the numbers 
and delay times which are associated with main events occurring before time 
t; however, when delays can be negative (or advances), we have also to consider 
delays associated with main events after time t. 

The behaviour at time t yields to an analysis of the pairwise behav- 
iour of (Y ,Y, ) for each main event in (0,t], There are four possibilities 
a b 

concerning the type~a and type-b events which terminate Y and Y, : both 

a b 

occur after time t, type-a occurs after time t and type-b before time t, 

vice-versa, and both occur before time t. Let there be n main events in 

(0,t], with and n - Z^ - Z^^^ - Z^^^ main events giving pairs 

of type-a and type-b events of the above four types, respectively. Let 

{Y (i),Y- (i)} for i = 1,2,...,Z be the delay times past t associated 
a b t 

(^ 

with main events counted by Z , and similarly let Y'(i) for i = 1,2 Z' 

tat 

and Y’(i) for i = 1,2,...,Z^^^ be the delay times past t asociated with 
b t 

the other two situations. 

The key feature of the present model which makes it tractable is that, 
given the n main events in (0,t], they are independent and uniformly dis- 
tributed over this interval. Thus we can see that 



P 



Z = k, Y (i) ^ a., Y^(i) b^ for i = 1,2, 
t a X b 1 

= k^, Y^(i) ^ a' for i = l,2,...,k^ 
Z^^^ = k^, Y^(i) ^ b’ for i = l,2,...,kj^ 



»k 



-yt (yt) 



n 



n! 



n=k+k 






n! (n-k-k^-k^)I k! k^! k^! 
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{t 1 Y^<v)dv} 



n-k-k^-k^ 



v=0 



l!jt j \"'^i)'3''} 



v=0 
t 



I V ''* '*''} 

V=0 

, t 

I P(Y^<v, Yj,>Wbpdv} . 



(3.1) 



In (3.1) the variable v Is used to denote that the typical main event being 
considered is located at time t - v, or a time point v before t. Letting 
t in (3.1) we obtain after a little simplification, and dropping the 

suffix t, the arbitrary time initial conditions as 



P 



Z = k,Y (i) a , Y (i) 2: b for i = 1,2, 

3. X D X 

= k^, Y^(i) ^ a’ for i = 1,2,.. .,k^ 
_^(b) = k^, Y^(i) ^ b’ for i = l,2,...,k^ 



,k 



= exp{- y J [l-P(Y^<v, Y^<v)]dv} 
v=0 



X (k!) ^ 71 [y 

k=l 



v=0 



-1 ^ 

X (k !)"'^ 77^ [y 

^ 1=1 



v=0 



-1 

X (k, !) ^ 7r*'[y 

^ i=l 



v=0 



P(Y >v+a^, Y,>v+b,)dv] 
a i b i 



P(Y^>v+a^, Y^<v)dv] 



P(Y <v, Y.>v+b!)dv]. 
a ’ b i 



(3.2) 
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Both (3.1) and (3.2) are valid for k,k^,k^ ^ 0 subject to the interpretation 
that when any of k,k^,k^ are zero the corresponding products are omitted. 

We next consider the behaviour at any particular time, which we take as 
zero, which is due to negative delays. We must consider the behaviour at 0 
due to main events in (0,t] and then let t The situation is similar 

to that just considered except that we are interested in the ^overshoots' of 
the delays past zero in the negative direction. For each event in (0,t] we 
either have a double overshoot, an overshoot of type-a only, an overshoot of 
type-b only, or no overshoot at all. Let these overshoots, measured as posi- 
tive quantities, be denoted by X^(i), Xj^(i) for i = l,2,...,k, 
i = l,2,...,k^ and X^(i) for i = l,2,...,k^. Here we have changed from 
Y to X in our notation of the 'forward' initial conditions, and will use 
5 instead of Z for counting the main events. The 'backward' initial condi- 
tions can now be obtained analogously to the forward ones, and yield after some 
reduction 

= k, x^(i) ^ a^, ^ for i = 1,2, ...,k 

= k , X'(i) S: a! for i = 1,2,. ..,k 
a’ a 1 ’ ’ ’ a 



(b) _ 



= k^, X^(i) ^ b^ for i = 1,2, ...,k^ 
0 



= exp{- y 



[1-P(Y >v, Y,>v)]dv} 
a D 



V = — CD 



X (k!) 



-1 



k 

TI [y 

i=l 



P (Y <v-a . , Y <v-b . ) dv ] 
a l b 1 



-1 ^ 

X (k !) "^ TT^ [y 
i=l 



X (k^!) ^ ^ [y 



0 



P(Y >v-a!, Y,_<v)dv] 
a lb 



v=— 0 ® 
0 



P(Y >v, Y, <v-b!)dv] . 
a b 1 



(3.3) 



V=-o 
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A point to notice about both the forward and backward initial conditions is 
that they do not simplify appreciably when the delays within pairs are inde- 
pendent. Even if and are independent in the model, the initial 

pairwise delays, there are any, will be dependent. 

4. The Numbers of Events in the Course of Being Delayed. 

One aspect of the bivariate delayed Poisson process, considered at an 

arbitrary time t, is the number of events in the main process generated 

before time t which are delayed until after time t, and conversely the 

number of events generated after time t which are advanced to before time 

t. Now the bivariate delayed Poisson process can be regarded as a simple 

bivariate case of Lewis’s branching Poisson process (Lewis, 1964), where a 

main event generates two subsidiary processes each of one event. In this 

terminology we will be obtaining the joint distribution of the numbers of each 

type of subsidiary process running at an arbitrary time. It is shown in Lewis 

(1964) for positive independent delays that marginally the distribution will 

be Poisson with means E(Y ), E(Y, ); and also, for positive independent 

a b 

delays, the joint Poisson distribution is stated in Cox and Lewis (1972). 

The result of two-sided dependent delays is derived here, and perhaps places 
the earlier results in a wider perspective. 

The required joint distribution is obtained from the forward and back- 
ward initial conditions by first setting the ^i’^i (3«2) and 

(3.3) equal zero. The backward and forward contributions will be independent 
because of the main Poisson process, and we require the joint distribution of 
{Z+Z^^\ Let the probability generating functions 

(p.g.f.’s) of the two pairs be and (p^(z^yZ^) . From (3.2) we thus 

have 
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oo 



OO 00 00 



<Pl(Zl»Z2) = I I I exp{- 



k=0 k =0 k, =0 

a v=0 



[1-P(Y <v, Y,<v)]dv 
a b 



X (k!) ^ {y I P(Y^>v, Y^>v)dv}^ 



v=0 



X (kgO 



P(Y >v, Y^<v)dv} 
a b 



v=0 



X (k^!) ^{y 



P(Y <v, Y >v)dv} 

H D 






v=0 



OO 



oo 



= exp{y(z^-l) 



P(Y >v)dv 4- (z -l)y 
J a 2 

v=0 



P(Y^>v)dv 

J b 

v=0 



+ y(z^-l) (z^-l) 



oo 

' 

P(Y >v, Y, >v)dv} . 

J a b 

v=0 



(4.1) 



This is the p.g.f. of a bivariate Poisson distribution X and is essentially 
the result stated by Cox and Lewis (1972, equation (4.34)). For the negative 
delays (or advances) we have from (3.3) 



= exp{y(z^-l) 



P(Y <v)dv + y(z„-l) 

Si z. 



P(Y^<v)dv 

b 



V=— oo 



Y =— 00 



4* y(zj^-l) (z^-l) 



P(Y <v, Y, <v)dv} . (4,2) 

a b 



Y=— oo 



The product of (4.1) and (4.2) gives the required p.g.f., and 

remains of the well known infinitely divisible bivariate Poisson form. The 

sum of the first integrals in (4.1) and (4.2) shows that the marginal means 

are E{|y |} and E{ |Y- |}, respectively. We now show that the covariance 
a D 

term 
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oo 



y{ P(Y <v, Y <v)dv + P(Y >v, Y, >v)dv}. 
Jab jab 

v=-a> v=0 



(4.3) 



is, apart from the factor y, the expectation of M(Y^,Y^) where 






for Y Y. s 0 
a b 



Min{|Y L |y^|} for Y Y, > 0. 
' a ' ’ ' b a b 



(4.4) 



This may be seen from the successive equalities 



E{M(Y ,Y, )} = E{min(Y ,Y,)|y >0, Y.>0}P(Y >0,Y.>0) 



a’ b' 'a 



a ’ b 



+ E{min(-Y ,-Y. )1 y <0, Y.<0}P(Y <0, Y,<0) 
a b a b a ’ b 



P(Y >v, Y >v Y >0, Y, >0)dv x P(Y >0, Y, >0) 
a’n'a b a b 



v=0 

’ 

j 

v=0 

* 

j 

v=0 



P(-Y >v, -Y^>v Y <0, Y^<0)dv x p(y <0, Y^<0) 
a’ b'a’b a’b 



P(Y >v, Y, >v)dv + 
a b 



P(Y <v, Y^<v)dv. 
a b 



(4.5) 



The interpretation of M(Y ,Y, ) is that it is the absolute length of the over- 

a b 

lap of a pair (Y ,Y, ); if Y and Y, are of opposite sign there is no 
a b a b 

overlap and M(Y^,Y^) is zero. The joint distribution of the numbers of 
events which are in the course of being delayed at an arbitrary time thus has 
joint p.g.f. 

(pCZj^fZ^) = exp{(z^-l)yE( [Y^|) + (z 2 -l)yE( [Y^ | ) + (z^^-l) (z 2 -l)yE[M(Y^,Y^) ] }. (4.6) 

When delays are non-negative, |Y^| = 



(4.6) reduces the result stated by Cox and Lewis (1972). We may also note 
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that (4.6) applies to a bivariate generalized branching Poisson process, 

Lewis (1969), if and are understood as the durations of subsidiary 

processes. Further to Cox and Lewis's remarks in section 4.6 of their 
Berkeley paper, it is interesting to note that the initial conditions cannot 
be constructed from the joint probabilities in (4.6) and the marginal forward 
recurrence time distributions of Y^ and Y^; this is because the dependency 
structure requires the joint distribution of these quantities. 



5. The Stationary Bivariate Counting Process. 

In this section we obtain the joint p.g.f. of 

the number of events of type-a and type-b which occur in the intervals (0,t ] 

a 

and (0,t^] of the stationary process. We know marginally that N^(t^) and 

N|^(t^) are Poisson random variables of means pt^ and so interest 

concerns their joint behaviour. Cox and Lewis (1972, Sections A. 3 and 4. A) 

obtained the joint distribution of N (t ) and N, (t, ) for pairwise independ- 

a a D D 

ent and non-negative delays, and we present here the corresponding result for 
pairwise dependent and two-sided delays; this more general case is shown to 
produce a rather simpler result. 

The type-a and type-b events in (0,t] arise from three sources: (1) 

those from main events in (-°°,0] which have been delayed into (0,t ] or 

a 

(0,t- ]; (2) those from main events in (0,max(t ,t, )] which have been delayed 
D a D , 

or advanced and are still in (0,t ] or (0,t, ]; (3) those from main events 

a D 

in (max(t ,t, ),<»] which have been advanced into (0,t ] or (0,t, ]. The 
a D a D 

contributions from these three sources are independently distributed in virtue 
of the Poisson main process, and we shall let their bivariate p.g.f 's be denoted 
by <p(z^,Z 2 ;-®, 0 ), and <p(z^, Z 2 ; t^, t^,“) , respectively. 



9 



Derivations of all three are similar; particular cases of the first two appear 
in Cox and Lewis (1972), and so we limit ourselves to the derivation of 
which arises from negative delays. We first state for 

future reference. 



cp (z^, Z 2 ;-“, 0 ) = exp{y(Zj^-l) P(v<Y^<v+t^)dv + p(z 2 -l) J P(v<Y^<v+t^)dv 

v=0 v=0 



+ y (Zj^-1) (Z 2 -I) 



P(v<Y <v+t , v<Y,_<v+t,_)dv} 
a a D b 



(5.1) 



v=0 



and 



cp(z^,Z2;0,t^,tj^) 



max(ta,tj^) 



max 



exp{y(z^-l) 



P(-v<Y <t -v)dv + y(z^-l) 
a a z 



P(-v<Yj^<t^-v)dv 



v=0 



v=0 



max(ta,tb) 



+ y(z^-l)(z2-l) 



P(-v<Y <t -V, -v<Y, <t, -v)dv>. (5.2) 

a a b b 



v=0 



For u)(z, ,z_;t ,t, ,«>) consider events which are advanced from (max(t ,t, ),T] 

1 z a 0 a b 

into (0,t^] or (0,tb], and ultimately let X For a main event at 

time max(t^,tb) + v, its contribution to cp (z^^, Z 2 ; t^» t^ ,") is 

1 + (z^-l)P(-v-tQ<Y^<-v-tb) + (z2-l)P(-v-tQ<Yb<-v-t2) 

+ (Zj^-1) (z2-l)P(-v-tQ<‘55^<-v-t^, -v-tQ<Yb<-v-t2) (5.3) 

where t„ = max(t ,t,), t, = max(t,-t ,0), t- = max(t -t, ,0). 
u abl ba z ab 



Conditional on the number of events in (tQ,X], the terms (5.3) are independ- 
ently and uniformly distributed over (tQ,I]. Hence, following the argument 
in section 3 to obtain the initial conditions, we obtain 
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oo 



oo 



= exp{u(z^-l) 



P(-v-tQ<Y^<-v-t^)dv + yCz^-l) 



P (-V- t^<Y^<-v- 12 > dv 



v=0 



v=0 



w 

+ y (Zj^-1) (Z2“l) I P(-v-tQ<Y^<-v-t^,-v-t|j<Y^<-v-t2)dv} . (5.4) 



v=0 



We observe that (5.1), (5.2) and (5.4) are all of the infinitely divisible 
bivariate Poisson form, and thus it remains only to obtain expressions for the 
parameters in their product; this is the required bivariate p.g.f. Simple 
changes of variables in the integral coefficients of (z^^-l) show that their 
sum may be written 



+00 



P(v<Y <v+t )dv, 
a a 



(5.5) 



If we assume that Y has a density then this integral evaluates to yt 

a a 

straight away; generally we will make this assumption. However, if Y is 

a 

restricted to non-negative values this may easily be proved without the assump- 
tion of a density. We thus have that N (t ) and N, (t, ) have marginal 

a a b D 

Poisson distributions of means yt^ and Similarly, adding together the 

coefficients of (z^-l)(z 2 -l) from (5.1), (5.2) and (5.4) gives 



-hco 



P(v<Y <v+t , 
a a 



v<Yb<v+t^)dv 



Y=— oo 



(5.6) 



as Cov{N^(t^), a quantity which is always non-negative. 

Thus, the general result for the bivariate p.g.f. of N (t ) and 






is 



cp(z^,Z 2 ;t^,t^) = exp{(z^-l)yt^ + (z 2 -l)yt^ 



+ (Zj^-1) (z2“l)y 



P(v<Y <v+t , v<Y <v+t, )dv} . (5.7) 

a a b b 



Y=— cx> 
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When Y and Y are non-negative, the covariance becomes 

3. D 



m: 




0<Y <t -V 


0 


0 


v<Y <v4-t 


y 


P 


a a 


dv + y 


p 


a a 


J 

v= 


=0 


0<Y, <t, -V 
D b 


« 

v= 


=0 





(5.8) 



and this makes (5.7) correspond to Cox and Lewis (1972, equation (4.25)). 
Setting = 22 = 0 in (5.7) gives the bivariate exponential forward 
recurrence times; further details are given by Cox and Lewis. 



6. The Semi-Synchronous Interval Distributions. 

Finally, in this paper we consider the distributions of intervals 

between events in the bivariate delayed Poisson process; the intervals between 

events of the same type have trivially, exponential distributions of parameter 

y. Consequently, we shall be interested in the distribution of intervals 

between consecutive events of opposite type; these have been called by Cox 

and Lewis (1972) the semi-synchronous intervals . We take an arbitrary event 

of type-a, say, and require the distribution of the interval to the first 

subsequent event of type-b; this interval will be denoted by Some 

a 

general theory relating synchronous, semi-synchronous and asynchronous interval 
distributions is given in Wisniewski (1972). We shall, however, develop our 
result from the definition in Cox and Lewis (1972), 



P(X^^^>x) 

a 



lira P{R (t,t+x) = 0|N (0,t) 1} , 

t+O ^ 



( 6 . 1 ) 



which is analogous to the definition in Khintchine (1960, Section 9) of Palm 
functions for synchronous interval distributions. Taking (6.1) we write it as 



P(X^*^^>x) = lim 
^ t+0 



P{Nj^(t,t+x) = 0} - P{N^(0,t) = 0, Nj^(t,t+x) = 0}' 



P{N (0,t) ^ 1} 
a 



( 6 . 2 ) 
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Dividing both the numerator and denominator by t, and assuming that the 
individual limits exist, as will be evident in a moment for the present 
process, we have the general result 



P(X^^^>x) = -y Vp{N^(0,t) = 0, N^(t,t+x) = (6.3) 

Here denotes a right hand derivative and y is the value of 

t ^P{N (0,t) ^1} as t i 0; the latter is trivially true from Khintchine’s 
a 

work since N (0,t) counts events in a Poisson process of rate y. 
a 

Our main concern now is thus to evaluate P{(N^(0,t) = 0, Nj^(t,t+x) = O) 

for the bivariate delayed Poisson process. Since N^(0,t) and Nj^(t,tH-x) 

refer to the stationary process we must consider main events which have occurred 

in the intervals (-T^,0], (0,t] and (t,t+T 2 ], eventually letting and 

A main event which occurs at -v in (-T^,0] does not lead to an a- 

event in (0,t] or a b-event in (t,t+x] when its delays are such that 

Y € [v,t+v]’, Y € [ t-Hv, t+v+x] ’ , where prime denotes the complementary 
a D 

interval; similarly for a main event at v in (0,t] the condition on the 

delays is Y € [-v,t-v]’, Y, € [ t-v, t-v+x] ’ , and finally for an event at 
a b 

t + V in (t,t+T^] the condition is Y € [- (t+v) ,-v] ’ , Y € [-v,x-v]’. 

an 

Use of the conditional property of the main Poisson process for each of the 
intervals (-T^,0], (0,t] and (t,t+T 2 ], leads to the result 



P{N (0,t) = 0, N, (t,t+x) = 0} = exp{y 
a D 



[-1+P 



Y €[v,t+v]' 'I 



Y €[t+v,t+x] ' 



] dv 



v=0 



L 

r 


rY €[-v,t-v]' ^ 


f 


rY 6[-(t+v),-v]''| 


J [-1+P 


a 

Y €[t-v, t-v+x] 'j 

D 


]dv + y [-1+P 


a 

Y^€[-v, x-v]' 



]dv} . (6.4) 



v=0 
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Slight changes of variable in the second and third terms then give 





OO 


rY ^[u,t+u]' 'i 


P{N (0,t) = 0, N, (t,t+x) = 0} = exp{y 
a D 


[-1+P 


a 

Yj^€ [t+u, t+u+x]j 



]du} . (6.5) 



U=— oo 



When the delays Y and Y, are non-negative, the corresponding result is 
a D 



P{N (0,t) = 0, N, (t,t+x) = 0 } = exp{y 
a b 



» 


fY €[u,t+u]’ ) 


[-1+P 

J 


a 

,Y € [t+u, t+u+x] 



u=0 



+ y 





rY >u ^ 


f 


[-1+P 

j 


a 

Y^6[u,u+x]\ 


]du - y 



]du 



P(Y,<u)du) • (6.6) 

b 



u=0 u=0 

We can now use (6.3) for the semi-synchronous interval distribution and have 



+00 



P{X^'^^>x} 



= 1 H 



U=-oo 



Y €[u,t+u]* 
a 

Y, €[ufu,t+u+x] 

D 









X exp{y 



[-1+P 



u= 



Y^€[u,t+u]’ >1 

Yj^€ [ t+u, t+u+x ] ' J ^ t=0 ) 



(6.7) 



The second term in (6.7) becomes exp(-yx) when t = 0 since the exponent 
is then of the form (5.5); the first term may be alternatively written as 



[P(u<Y <t+u) + P(t+u<Y,< t+u+x) - P 
a b 



u=— 00 



fu<Y <t+u ^ 
b 



a 

t+u<Y, 



<fh.+K]‘‘“il<t-o) ’ 



which then simplifies to 





r +<» 


+ 






1 - D 

t 




P 




ki= 


=—00 



u<Y <t+u 
a 



t+u<Y^ < t+u+xj 



du 



(t=0) • 



(6.8) 
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After writing the derivative in the incremental way and examining the expression 

at t = 0, it is found that (6.8) becomes P(Y -Y [0,x]*). Thus from (6.7), 

b a 

we have the main result of this section* as 



P{X^^^>x} = P(Yj^-Y^€[0,x]')exp(-yx), 



(6.9) 



and it is clear that this is a proper survivor function. One particular case 
of (6.9) is when Y^ and Y^ have exponential distributions of parameters 
and X^; then the p.d.f. of the mixed-exponential form 



,(b) , , \ 

A +X. A +A, 



-yx 



( 6 . 10 ) 



a b 



a b 



The distribution of Y, - Y appearing in (6.9) also features in other 

b a 

aspects of the bivariate delayed Poisson process, as mentioned by Cox and Lewis 
(1972), and the reason is fairly evident. For some purposes, the main Poisson 
process can be eliminated from the model by concentrating on (say) just the 
type-a events; a type-b event occurs at a distance Y^ - Y^ from its type-a 
event, and these follow a Poisson process anyway. Such reasoning leads to a 
quicker derivation of (6.9), but our aim here has been to apply, for the first 
time it is believed, the semi-synchronous definition (6.1). It should be noted 
however, that the univariate process with displacements Y^ - Y^ is not iden- 
tical to our delayed Poisson process in all respects; for instance, in terms 
of initial conditions and numbers of events in the course of being delayed. 

Some joint counting distributions for the univariate delayed Poisson process 
have been obtained by Milne (1970) when pursuing identif iability questions. 

r- • • \ r \r (^ ) 

Finally, we note that a joint distribution of interest is that of 
(a ) 

together with that of , the interval from the arbitrary type-a event to 

3 . 
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the next type-a event. We can define generally the required joint distribution 
as 

P{X^®^>x, X^^^>y} = lim P{N (t,t+x) = 0, N^(t,t+y) = 0 | N (0,t) 2: 1> , (6.11) 

a a I ^ ^ D a 

t'l^O 



or use the univariate delayed Poisson process with ^ ^ ” ^a dis- 

placements. We follow the latter course. Taking our origin as an event in the 
Poisson process of type-a events, we require that there be no type-a events in 
(0,x] and no type-b events in (0,y]. We consider type-a events which have 
been generated in (-Tj^,0] and (x,!^] and then let and T 2 tend to 

infinity; thus 



P(X^^^>x, X^^^>y) = 
a a 



lim I 
r=0 

T -»<» 

2 



® -ViT. (yT. r 



1 1^ 



r ! 



-^1 



P(Y6 [v, vfy ] * )dv x e x 






-Px 



v=0 



“ -PT (PT ) r, T 2 -IS 

I e — -J — P(Y€[-v-x,y-v-x] ')dv P(Y€[0,y]') 

s=0 '-^2 

v=0 



(6.12) 



which may be simplified to 

0 

P(X^^^>x X^^^>y) = exp{-yx-yy+y [ P(u<Y -Y <u-Hy)du} x P (Y -Y €[0,y]^). (6.13) 

a a I D a Da 

u=-x 

This joint distribution has the marginal distributions already discussed. From 
a statistical analysis point of view it would be useful to obtain the correla- 
tion between X^^^ and X^^^; however, it appears from (6.13) that no simple 
a a 

general expression will be available, although (6.13) should yield the result 
with particular assumptions on the delay structure. 
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